Ingredients of a Casimir Analog Computer 
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We present the basic ingredients of a technique to compute quantum Casimir forces at micrometer 
scales using antenna measurements at tabletop (e.g. centimeter) scales, forming a type of analog 
computer for the Casimir force. This technique relies on a correspondence that we derive between 
the contour integration of the Casimir force in the complex frequency plane and the electromagnetic 
response of a physical dissipative medium in a finite real-frequency bandwidth. 
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Casimir forces arise due to quantum fluctuations of the 
electromagnetic field pQ and can play a significant role in 
the physics of neutral, macroscopic bodies at micrometer 
separations, such as in new generations of microelectronic 
mechanical systems (MEMS) |2j [3]. These forces have 
previously been studied both in delicate experiments at 
micron and sub-micron lengthscales [4 and also in theo- 
retical calculations that are only recently becoming fea- 
sible for complex non-planar geometries [3[6]. Here, we 
propose a third alternative by deriving an equivalence be- 
tween quantum fluctuations and the classical electromag- 
netic response in bodies separated by a conducting fluid, 
as illustrated in Fig. [I] Using this equivalence, we pro- 
pose the possibility of experimental /S-matrix measure- 
ments for microwave antennas in centimeter-scale models 
that indirectly yield the Casimir force between micron- 
scale objects. Such a centimeter-scale model is not a 
Casimir "simulator," in that one is not measuring forces, 
but rather a quantity that is mathematical related to 
the micron-scale Casimir force — in this sense, it is a kind 
of analog computer. We believe that this mathematical 
equivalence between disparate quantum and classical sys- 
tems reveals new opportunities for the experimental and 
theoretical study of Casimir interactions. 

In the following letter, we first review a well-known 
formulation of the Casimir force in terms of the classi- 
cal Green's function (GF) via the electromagnetic stress 
tensor (ST) and the fluctuation-dissipation theorem [7]. 
Although this formulation is normally expressed for ei- 
ther real or imaginary frequencies, we consider the gen- 
eral complex-frequency (lj) plane. We then show that 
the mapping to complex frequency is equivalent to a 
real-frequency GF with a transformed electromagnetic 
medium e c . From this point of view, however, it turns 
out that only certain contours in the complex-frequency 
plane correspond to physically realizable e c , and in par- 
ticular the real and imaginary frequency axes are un- 
suitable. Instead, we identify another contour and show 
its equivalence to a conventional conducting medium, 
demonstrate that the response of such a medium yields 
the correct Casimir force in nontrivial geometries, and 
consider the implications and possible materials for prac- 
tical experiments. The key point is that, once the 
Casimir force is expressed in terms of the response of a 
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FIG. 1: Schematic illustration of correspondence between ob- 
servations of Casimir forces at fim. lengthscales and antenna 
measurements of the electromagnetic response functions at 
tabletop lengthscales (e.g. cm). 



realizable medium over a reasonably narrow bandwidth, 
the scale-invariance of Maxwell's equations permits this 
response to be measured at any desired lengthscale, e.g. 
in a tabletop microwave experiment. 

The Casimir force can be expressed as an integral of 
the mean electromagnetic ST over all frequencies [7]. The 
mean ST is determined simply from the classical GF (the 
fields in response to current sources at a fixed frequency) , 
thanks to the fluctuation-dissipation theorem. It turns 
out, however, that this frequency integral is badly be- 
haved from the perspective of numerical calculations (or 
experiments, below). Fortunately, because the integrand 
is analytic, one can deform the integration contour into 
the complex-frequency uo plane. 

More generally, given an arbitrary contour (for 
convenience below, we parameterize the contour by a real 
£), the force in the z-th coordinate direction is given by: 
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The standard Wick rotation corresponds to the partic- 
ular choice uj{£) = i£ and yields a smooth and rapidly 
decaying integrand [5 . The mean ST (T^) is related to 
the electric (E) and magnetic (H) field correlation func- 
tions by the standard equation (assuming non-magnetic 
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materials, /i = 1, for simplicity): 
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The field correlation functions are, in turn, related to the 
frequency-domain classical photon GF, Gij(uo; r, r'), by 
the fluctuation-dissipation theorem: 

(E^E^v')) = -u 2 G ij {u-v,v') (3) 

7T 

(Hi^Hjir')) = --(Vx)«(V'x) jm G <TO (a;;r,r'), (4) 
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where Gij satisfies Maxwell's equations: 

[V x V x - e(r, uo)uo 2 ] Gj (uo; r, r') = tf (r - r')^ (5) 

Equation ([5| can be solved in a number of ways, for ex- 
ample by a finite-difference discretization [5] or even an- 
alytically in one dimension [7 j. Of course, the diagonal 
(V = r) part of the GF is formally infinite, but this 
singularity is not relevant because its surface integral is 
zero, and it is typically removed by some regularization 
(e.g. by the finite discretization or by a finite antenna 
size in the proposed experiments below). A crucial step, 
as mentioned above, is the passage to imaginary frequen- 
cies — i£. For real frequencies, the GF is oscillatory, 
leading to a highly oscillatory ST integrand that does not 
decay — even when a regularization (ultraviolet cutoff) is 
imposed, integrating a highly oscillatory function over a 
broad bandwidth is problematic. For imaginary frequen- 
cies, on the other hand, the GF is exponentially decaying, 
due to the operator in Eq. ([5| becoming positive-definite 
(V x V x + s£ ) 2 ) [5], leading to a decaying non-oscillatory 
integrand. 

However, the Wick rotation is not the only contour 
in the complex plane that leads to a well-behaved de- 
caying integrand. This is illustrated by Fig. [2j which 
shows the force integrand in the complex plane for 
the piston-like geometry of Fig. [I] In particular, we 
calculate the (x-direction) force integrand dF x /duo = 
<£j> surf T x j(r, uj)dSj on one square, for geometric pa- 
rameters h = 0.5d and s = where d is the separa- 
tion between the blocks. Here, we plot ln\RedF x /duo\, 
which illustrates the basic features of the integrand 
duo / d£ dF x / duo (lmdF x /duo is also important, but is qual- 
itatively similar). As described above, this integrand is 
oscillating along the Re a; axis and decaying along the 
Imuo axis, but it also decays along any contour where 
Imuo is increasing (such as the three contours shown, to 
be considered in more detail below). 

The fact that uo and e only appear together in Eq. ([5]), 
as euo 2 , immediately suggests that, instead of changing 
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FIG. 2: Complex-frequency uj plot of the Casimir force inte- 
grand (In | Re dF x /duo \ ) , where dF x /dcu is in units of h/d 2 , for 
the geometry of Fig. [l] As the real-uu axis is approached, the 
integrand becomes highly oscillatory, which is only partially 
revealed here due to the finite frequency resolution. Vari- 
ous integration contours of interest are labelled as black and 
dashed lines. (Inset:) Vacuum e — 1 contour deformations 
and their corresponding (real- frequency) physical real- 
izations £ c (0- 



uo to a complex number, we can instead operate at real 
frequencies by transforming e. In particular, a complex 
is equivalent to a real frequency £ with a complex 
£, where the imaginary part of e corresponds to dissipa- 
tion loss. Thus, an intuitive explanation for why trans- 
forming to the complex-frequency plane was numerically 
useful above is simply that it corresponds to lossy ma- 
terials that damp out the oscillations. Given a medium 
e at a complex frequency uo{£) where we wish to com- 
pute the ST for the Casimir force, we will transform to 
an equivalent problem at a real frequency £ and com- 
plex permittivity e c . Namely, operating at a complex 
frequency a;(£) is clearly equivalent (for the photon GF) 
to operating at a real frequency £ and transforming a 
given material via — > e(uo)uo 2 /£ ) 2 = e c . Conversely, 
any frequency-dependent material £ c (0 a t a given point 
in space can be related to the GF for vacuum (e = 1) 
at that point by going from the real frequency £ to a 
complex frequency uo = £y/e c (£)- 

Because the ST is expressed in terms of the GF, and 
the GF at microwave lengthscales is merely a rescaling 
of the GF at micron lengthscales (if suitable materials 
can be found), one can conceivably measure the GF in 
an experiment via the ^-matrix elements of antennas at 
centimeter scales, and so determine the Casimir force via 
integration of the ST. The passage to complex frequencies 
is essential here, as in numerics, because measuring the 
real-frequency GF will yield a highly oscillatory force in- 
tegrand over an infinite bandwidth, imposing significant 
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experimental challenges. Unfortunately, it is difficult 
to implement complex frequency deformations directly, 
because a complex frequency corresponds to fields and 
sources with exponential growth in time. An alternative, 
suggested by the correspondence above, is to measure the 
real-frequency GF of a physical medium with a complex 
permittivity £ c (£) = e{uo)uo 2 / '£ 2 . This e c medium should 
satisfy two properties: it should correspond to an uo con- 
tour where the ST integrand is rapidly decaying, and it 
should be physically realizable. For simplicity, we begin 
by considering complex contours for the ST in vacuum 
(e = 1). The extension to arbitrary geometries/materials 
is straight-forward: an arbitrary inhomogeneous medium 
e(r, £) corresponds to e c = uu 2 e(r, oo)/^ 2 . 

If £ c (£) is to correspond to a physical medium, it must 
satisfy the complex-conjugate property £ c (— £) = £ c (0* 
as well as the Kramers-Kronig (K-K) relations [8 j. It is 
most important to satisfy these conditions for small £, 
since the ST integrand is dominated by long- wavelength 
contributions. One should also prohibit gain media, 
which would lead to the exponentially growing fields we 
are trying to avoid by not using complex uo. For ex- 
ample, Wick rotations correspond to £ c (£) = —1, and 
this is only possible at £ = in a gain medium, since 
in a dissipative medium, e c is real and positive along 
the whole imaginary-^ axis (this is implied by K-K). 
The generalization to arbitrary rotations in the com- 
plex plane e = e 2 ^ is both a gain medium and violates 
£ c (— £) = £ c (0* near £ = 0- Thus, no realizable material 
can emulate these contours even in a narrow bandwidth 
around £ = 0, as summarized on the table (inset) of 
Fig-H 

Although traditional Wick rotations correspond to un- 
physical materials, there are obviously many physical 
lossy materials to choose from, each of which corresponds 
to a contour in the complex plane, and one merely needs 
to find such a "physical" contour on which the ST is 
rapidly decaying so that experiments can be performed 
over reasonable bandwidths. A simple and effective lossy 
material for this purpose is a conductor with conductiv- 
ity a. Because the integral will turn out to be dominated 
by the contributions near zero frequency, it is sufficient 
to consider a to be a constant (the DC conductivity), al- 
though of course the full experimental permittivity £ c (£) 
could also be used. Specifically, consider the general class 
of conductors defined by dispersion relations of the form 
£ c (0 = 1 + 2cr/£, corresponding to vacuum with a com- 
plex contour uo a = ^s/i + icr/£. As shown in Fig. |5J the 
integrand of this contour is in fact well behaved, rapidly 
decaying and exhibits few oscillations. 

We now consider the Casimir force for the same struc- 
ture as in Fig. [TJ still calculated by the same finite- 
difference method as in Fig. [3j but we now focus on the 
properties along different contour choices for both physi- 
cal and unphysical media. In particular, Fig.[2|top) plots 
the partial integral {duo / d£) (dF x /dw)d£, normalized by 
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FIG. 3: (Top:) Partial force integral J Q dF x , normalized by 
F x , as a function of £, for the various cj(£)-contours (equiv- 
alently, various e c — uj 2 /^ 2 ) shown in Fig. [3] The solid 
green, red and black lines correspond to conductive media 
with a = 10, 10 2 and 10 3 , respectively (a has units of c/d). 
The dashed grey and solid blue lines correspond to <j) — tt/4 
and <fr = tv/2 (Wick) rotations. (Bottom:) Illustration of the 
required frequency bandwitdth for a possible realizations us- 
ing saline solution at separation d = 30cm. The red lines plot 
the xx-component of the photon GF G xx at a single location 
on the surface contour (see inset) as a function of £ (GHz). 
The black line is the corresponding partial force integrand. 

the total force J °° dF x , as a function of £. [As it must, 
the total integral over £, the force F x , is invariant re- 
gardless of the contour a;(£) and agrees with previous 
results [5 ; Specifically, F x = 0.0335 (hc/d 3 ).} We now 
comment on two important features of the uo a contour 
that are relevant to experiments. 

First, the Jacobian factor for uo a is given by duo a /d^ = 
0.5(2 + icr/£) j y/l + icr/t; and turns out to be very im- 
portant at low £. The ST integrand itself goes to a con- 
stant as £ — > (due to the constant contribution of zero- 
frequency modes) , but the Jacobian factor diverges in an 
integrable square-root singularity ~ VW£- Since this 
singularity is known analytically, however, separate from 
the measured or calculated GF, integrating it accurately 
poses no challenge. Second, the larger the value of a, the 
more rapidly the ST integrand decays with £, and as a 
consequence the force integral for larger a is dominated 
by smaller £ contributions. In comparison, previous cal- 
culations of Casimir forces along the imaginary-frequency 
axis revealed that the relevant £ bandwidth was deter- 
mined by some characteristic lengthscale of the geometry 
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such as body separations [5 j. Here, we have introduced 
a new parameter a that can squeeze the relevant £ band- 
width into a narrower region. This "spectral squeezing" 
effect is potentially useful for experiments, as it partially 
decouples the experimental lengthscale of the geometry 
from the required frequency bandwidth. 

As a consequence of the above results, we can now 
outline a possible experiment at centimeter lengthscales 
that determines the Casimir force at micron lengthscales, 
a Casimir analog computer (CAC). Suppose that one 
wishes to compute the Casimir force between perfect- 
metal objects separated by vacuum, such as the geome- 
try in Fig. [I] One would then construct a scale model 
of this geometry at a tabletop scale (e.g., centimeters) 
out of metallic objects (which can be treated as perfect 
metals at microwave and longer wavelengths). To deter- 
mine the ST integrand along a complex-ccv contour, one 
would measure the GF at real frequencies £ for the model 
immersed in a conducting fluid. The GF is related to the 
S- matrix of pairs of antennas, and the diagonal of the GF 
to the ^-matrix diagonal of a single antenna [noting that 
the finite size of the antenna automatically regularizes 
the integrand, as noted after Eq. It is important for 
the model structure to be large enough that the introduc- 
tion of a small dipole-like antenna does not significantly 
alter the electromagnetic response. In general, the ST 
must be integrated in space over a closed surface around 
the object, and correspondingly the antenna's ^-matrix 
spectrum must be measured at a number of antenna po- 
sitions (2d quadrature points) around this surface. (Un- 
less one is interested in computing the force on a sin- 
gle atom, which requires a single antenna measurement.) 
The different components of the GF tensor correspond 
to different antenna orientations. The magnetic GF can 
be determined from the photon GF by Eq. ([3]), or possi- 
bly by employing "magnetic dipole" antennas formed by 
small current loops. 

We now consider a particular CAC (at the cm scale) 
that employs realistic geometric and material parame- 
ters. Many available fluids exhibit almost exactly the 
desired material properties from above. One such exam- 
ple is saline water, which has = e s + ia/%, where 
e s ~ 80 and a ~ 5 S/m for relatively small values of 
salt concentration [9]. A calculation using these param- 
eters, based on the geometry of Fig. [I] assuming object 
sizes and separations at the centimeter to meter scale 
(we choose d = 0.3 m for the structure in Fig. [I] cor- 
responding to a frequency of 1 GHz), reveals that it is 
only required to integrate the stress tensor up to small 
GHz frequencies £, which is well within the reach of con- 
ventional antennas and electronics. This is illustrated 
in Fig.[3jbottom), which plots the G xx component of the 
GF (red lines) as well as the partial force integrand (black 
line), showing the high £ < 1 GHz cancellations that 
occur once the ST is integrated along a surface (inset). 
We note that most salts exhibit additional dispersion for 



£ > 10 GHz [9], but we do not need to reach those fre- 
quency scales. (Nevertheless, should there be substantial 
dispersion in the conducting fluid, one could easily take 
it into account as a different complex-o; contour.) 

Some attention to detail is required in applying this 
correspondence correctly. For instance, using network 
analyzers, what is measured in such an experiment is not 
the photon GF G, but rather the electric ^-matrix 
(the currents in a set of receiver antennas due to cur- 
rents in the source antennas) , related to the electric GF 
(the E-field response to an electric current J) by a fac- 
tor depending on the antenna geometry alone (relating 
J to E). The electric GF will differ from the photon GF 
by a factor of the real frequency To summarize, the 
photon GF will be given in terms of the measured by 
Gij(uj) = (a/i£)Sij(£), where a is the antenna-dependent 
geometric factor. To obtain the ST from G^-, one multi- 
plies by factors of as in Eq. H}. 

Figure ^bottom) 

illustrates the expected behavior of 

$xx ^ G X x I £ in & re- 
alistic system employing a sline solution with d = 30cm. 

The use of tabletop models and analog computers in 
physics, though rarely explored in the context of vac- 
uum fluctuations, continues to play an important role in 
contemporary research fields, such as quantum informa- 
tion [10 . Especially for three-dimensional geometries, 
tabletop experiments offer a route to rapidly exploring 
many different geometric configurations that remain ex- 
tremely challenging for conventional numerical calcula- 
tion. Although many details of such an experiment re- 
main to be developed, we believe that the basic ingredi- 
ents are both clear and feasible, at least when restricted 
to perfect-metal bodies. The most difficult case to realize 
seems to be the force between imperfect-metal or dielec- 
tric bodies with a permittivity e(u), as the corresponding 
tabletop system requires materials with a specified dis- 
persion relation e^ ody (f) / 'e^ uid (f) = e(u(£)) relative to 
the conducting fluid. This may be an opportunity for 
specially designed meta-materials with the desired fre- 
quency response. 
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